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Abstract. We revisit the Near Equidiffusional Flames (NEF) model intro- 
duced by Matkowsky and Sivashinsky in 1979 and consider a simplified, quasi- 
steady version of it. This simplification allows, near the planar front, an ex- 
plicit derivation of the front equation. The latter is a pseudodifferential fully 
nonlinear parabolic equation of the fourth-order. First, we study the (orbital) 
stability of the null solution. Second, introducing a parameter e, we rescale 
both the dependent and independent variables and prove rigourously the con- 
vergence to the solution of the Kuramoto-Sivashinsky equation as e — » 0. 



Flames constitute a complex physical system involving fluid dynamics, multistep 
chemical kinetics, as well as molecular and radiative transfer. The laminar flames 
of low-Lewis-number premixtures are known to display diffusive-thermal instability 
responsible for the formation of a non-steady cellular structure (see [H]). However, 
the cellular instability is quite robust against these aero-thermo-chemical complex- 
ities and may be successfully captured by a model involving only two equations: 
the heat equation for the system's temperature and the diffusion equation for the 
deficient reactant's concentration. In suitably chosen units, the so-called thermal- 
diffusional model reads, see e.g., [6]: 
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(1.3) 



(1.2) 



Yt = Le-\Y,, + Yyy) - niY, 8), 
n = ^Le-^P^Yexp[f3{e - l)/{a + (1 - (7)6)]. 



Here, 8 = (T — Tu)/{Tad — r„) is the scaled temperature, where T„ and Tad cor- 
respond to the temperature of the unburned gas, and the adiabatic temperature 
of combustion products, respectively; Y = C/Cu is the scaled concentration of the 
deficient reactant with Cu being its value in the unburned gas; x, y, t are the scaled 
spatiotemporal coordinates referred to Dth/U, and Dth/U'^, respectively, where Dth 
is the thermal diffusivity of the mixture and U is the velocity of the undisturbed 
planar flame; Le is the Lewis number (the ratio of thermal and molecular diffu- 
sivities); a = Tu/Tad] P = Ja(l — cr)/Tad is the Zeldovich number, assumed to be 
large, where Ta is the activation temperature; is the scaled reaction rate, where 
the normalizing factor ^Le^^0^ ensures that at /? ^ 1 the planar flame propagates 
at the velocity close to unity. 

Due to the distributed nature of the reaction rate il. Equations (|1.1[) and (|1.2[) are 
still difficult for a theoretical exploration. One therefore turns to the conventional 
high activation energy limit (/? ^ 1) which converts the reaction rate term into 
a localized source distributed over a certain interface x = S.(t,y), the flame front. 
Intensity of the source varies along the front as exp (|(8/ — 1)) (see [H]). Here, 8/ 
is the scaled temperature at the curved front, which may differ from unity (T = Tad) 
by a quantity of the order of f3~^. Due to the strong temperature dependence of the 
reaction rate {f3 ^ 1), even slight changes of 8/ may markedly affect its intensity, 
and thereby also local flame speed. The study of flame propagation is thus reduced 
to a free-interface problem. To ensure that the emerging free-interface model does 
not involve large parameters one should combine the limit of large activation energy 
(/3 3> 1) with the requirement that the product a = ^f3{l — Le) remains finite, i.e., 
the ratio of thermal and molecular diffusivities {Le) should be closed to unity. This 
is the Near Equidiffusive Flames model, in short NEF, introduced in [Kj. As a 
result, instead of the reaction diffusion problem for 8 and Y, one ends up with 
a free-interface problem for the new scaled temperature 9 = lim^^+oo 8 and the 
reduced enthalpy S — lim^^+oo (3~^{Q + Y — 1). More precisely, the system for the 
temperature 9, the enthalpy S and the moving ffame front, defined hy x = S,{t,y), 
reads 



For some mathematical results about this problem, see [5j [9l [TOl [TTl [121 [7] . Here, 
we consider only the case when a is positive, i.e., Le < 1. It will be convenient 
to assume periodicity in y with period £, and restrict attention to y € [—i/2,i/2]. 
At the front, 9 and S are continuous, the following jump conditions occur for the 
normal derivatives: 



(1.4) 



— = A0, x<at,y), 
9^1, x>at,y), 

— =AS-aA9, x^i{t,y). 



(1.5) 



(1.6) 



(1.7) 



exp(S'), 



dn 
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System (|1.4|) - (|1.8[) admits a planar travelling wave (TW) solution, with velocity — 1, 



d{x) 



expx, 
1, 



a: < 0, 
a; > 



Six) 



ax exp X, 
0, 



a; < 0, 
a: > 0. 



As usual one fixes the free boundary. We set ^{t, y) 
In this new framework: 



't + (p{t,y), x' = x-^{t,y). 



(1.9) 
(1.10) 
(1.11) 
where 



9t + ii~ ipt)e^' - A^e, x' < 0, 
e{x') = 1, x' > 0, 

St + il-^t)S.,, = A^S~aA^e, 



x'^0, 



A^ = (1 + {'■Pyf)Dx'x' + Dyy — iPyyDx' — 2LPyDx'y 

The front is now fixed at x' = 0. The first condition (|1.7p reads: 



the second one (jl.Sp becomes 

dS 



dx' 



= -exp(S'), 



de 



dx' 



A very challenging problem is the derivation of a single equation for the interface 
or moving front Lp, which may capture most of the dynamics and, as a conse- 
quence, yields a reduction of the effective dimensionality of the system. In this 
spirit, one of the authors in [17] derived asymptotically from the System (|1.4|) - (|1.8p 
the Kuramoto-Sivashinsky (K-S) equation in rescaled dependent and independent 
variables 

1 



(1.12) 



0. 



Since then, this equation has received considerable attention from the mathematical 
community. We refer to the book [20] and its extensive bibliography. 

This paper is devoted to a quasi-steady version of the NEF model. As a matter 
of fact, it has been observed in similar problems (see [^) that not far from the 
instability threshold the time derivatives in the temperature and enthalpy equations 
have a relatively small effect on the solution. The dynamics appears to be essentially 
driven by the front. Based on this observation one can define a quasi-steady NEF 
model replacing (|1.9[) - (|1.11[) by 

{l-ift)0^ A^e, a;'<0, 
e = l, a;' > 0, 

{l~ift)S.x. ^ A^S~aA^e, x'^0. 
Next we consider the perturbations of temperature u and enthalpy v. 

e^0 + u, S = S + v. 
Writing for simplicity x instead of x' , the problem for the triplet (u,v,ip) reads: 
(1 - ipt)ux - A^u - ipt6x = (Ay, - A)^, X < 0, 
M = 0, a: > 0, 
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(1 - Lpt)vx - ^^{v - au) - iptSx = [A^ - A){S - ad), x 7^ 0, 

where 

As in [31 13] , we introduce further simphfications: we keep only linear and second- 
order terms for the perturbation of the front (p, and first-order terms for the per- 
turbations of temperature u and enthalpy v. This leads to the equations: 

ux- Au~ iptOx = {A^ " A)e, a; < 0, 

Vx - A{v - au) - iptSx = (A^ - A){S - a9), x ^ 0. 
At a; = there are several conditions. First 

M = [v] = 0, 

however, since u{x) — for a; > 0, this is equivalent to 

u{0-) = [v] = 0. 

Second, 



'l + iipy)^0x + ux] = -exp(S' + w), 
hence up to the second-order: 



-1 + K] ^ -(1 + i^y)')-^ ^ 1 - -i^y)' ) ( 1 + t;(0) + -(z;(0))^ 



1 

and keeping only the first-order for v yields: 

[vx] = -aux{0). 
Therefore, the final system reads: 

Ux- Au - iptOx {{iPy)'^ - 'Pyy)dx, X < 0, 

Vx - A{v - au) - iptSx = {ipy)'^Sx - 'PyyS, x 7^ 0, 

(1.13) <j u{0) = [v]^0, 

viO) - Ux{0) = ^{iPy)^ 

[vx] = -aux{0). 

We remark that the equation for u associated with the boundary condition m(0) — 
entirely determines u when ip is given. Therefore, it can be viewed as a kind 
of pseudo-differential Stefan condition. We will take advantage of this remark in 
Section [31 

The goal of this paper is to show that this simplified NEF model still contains 
the dynamics of the system. It is simple enough to be integrated explicitly via a 
discrete Fourier transform in the variable y and therefore it allows a separation of 
the dependent variables. We get to a self-consistent pseudo-diff'erential equation for 
the front (p which reads: 

+ aXk - a)(pt{t, k) = (-4A2 + (a - l)Afe)<^(t, k) 

(1.14) +-^{Xl-iXl-AaXk + Aa)J^{t,k), A; = 0, 1, . . . , 
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where the — A^'s are the non-positive eigenvalues of the operator Dyy with periodic 
boundary conditions at y = zL£/2 (that we denote below by A) and 



is the symbol of operator — 4,Dyy. 

Equation (fTM]) can be written in the more abstract form: 

(1.15) |^ = J^(<^)+^((^^)2)^ 

where ^ is a pseudodifferential operator whose leading part is Dyy and G is a 
nonlinear operator whose leading term is j ^1 — 4:Dyy. This makes (|1.15p a strongly 
nonlinear equation, more precisely it is a fully nonlinear parabolic equation: in the 
L^-setting the nonlinear part is exactly of the same order as the linear operator. 
This is one of the main issues of this paper. Note that the realization of the operator 
y^T^^TD^ in the space of continuous and ^-periodic functions (say Cj) is defined 
only in a proper subspace of (the space of all the ^-periodic C-'^-functions) . Hence, 
in the Cj-setting, the nonlinear term J^(((^y)^) represents the leading part of the 
right-hand side of p.lSp . This would make the study of (jl.lSp more difhcult than 
in the L^-setting, where we confine our analysis. 

In the case when (p is smoother, we can rewrite Equation (jl.lSp as a fourth-order 
equation as follows: 

(1.16) ^^^-^(^) + ,^((^,)2), 

where is nothing but the usual fourth-order differential operator 

y{(p) = -tPyyyy ~ {a ~ l)'Pyy. 

Operators and ^ are pseudo-differential ones with symbols, respectively, 



bk ^Xi + aXk-a, fk = \ {Xl - iXl - AaX^ + 4a) . 



Therefore, 



^I-ADyy+aUl-ADyy-I 



^Dyy)i ADyy) ~ « ( ^ 1 " yy - 

The main feature of Equation (|1.16p is that the nonlinear part is rather unusual. 
Actually, it has a fourth-order leading term, as ,5^ has. Therefore (|1.16p is also a 
fully nonlinear problem, in contrast to (jl.l2p which is semilinear. 

The paper is organized as follows: In Section [3] we derive the front equation via 
an explicit computation of (|1.13p in the strip R x [—1/2,1/2]. Then, in Section [H 
we prove the following result 



Theorem 1.1. Let 

(1.17) a, = 1 



167r2 



Then, the following properties are met. 

(a) If a < ac, then, the null solution to Equation (jl.l6p is [orbitally) stable, with 
asymptotic phase, with respect to sufficiently smooth and small perturbations. 

(b) If a > ac, then the null solution to Equation (11.16^ is unstable. 
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An important question, that we address in Section [51 is the hnk between (|1.16p 
and K-S. FoUowing J/7j, we introduce a smaU parameter e > 0, setting 

« = ! + £, 

and define the rescaled dependent and independent variables accordingly: 

We see that solves the equation 

^ {/ - AeDrin + (1 + e) [^/I - ^eD^r, - l) } V' 

+ i{(/-4ei^^^)i -3(/-4ei^^^)-4(l + e) (^1 - 4eA,„ - /) } (A^V-)'- 

Then, we anticipate, in the limit e — > 0, that ijj ^ where $ solves (|1.12p . More 
precisely, we take for t. 

which blows up as e ^ 0; hence ac = 1 + ''4i- £■ Thus, ^Q becomes the new 
bifurcation parameter. We shall assume that i'o > 47r in order to have ac G (1, 
i.e., a > ac, otherwise the trivial solution is stable and the dynamics is trivial. 
The second main result of the paper is the following. 

Theorem 1.2. Let $0 e H"^ be a periodic function of period £q. Further, let $ 
be the periodic solution of (|1.12p (with period £q) on a fixed time interval [0,T], 
satisfying the initial condition $(0, •) — ^q. Then, if m is large enough, there exists 
£0 = £o{T) g (0, 1) such that, for < e < Eq, Problem (jl.l6p admits a unique 
classical solution ip on [0, T je^], which is periodic with period £qI \fs with respect to 
y, and satisfies 



V3(0,y) =£$o(2/%/i), |y| < 



4 



Moreover, there exists a positive constant C , independent of e € (0,eo]; such that 
\ip{t,y)-e<^ite^,yV^)\<Ce^, 0<t<^,\y\< 



e2' - 2^' 
for any e £ (0, Col- 
in other words, starting from the same configuration, the solution of p.l6p re- 
mains on a fixed time interval close to the solution of K-S up to some renormal- 
ization, uniformly in e sufficiently small. Note that the initial condition for ip is of 
special type, compatible with $0 and (|1.12p at t = 0. Initial conditions of this type 
have been already considered in [ij 12 12 • 

Although energy methods are known to be usually inefficient in fully nonlinear 
problems, here we may take advantage of the special structure of It allows us to 
establish sharp a priori estimates on the remainder (more precisely on its derivative) 
when £ is small enough. A key point is an extension of a lemma that we already 
successfully used in [U |4] . 

Finally, in the Appendix, for the reader's convenience we provide a quite detailed 
proof of the existence, uniqueness and regularity of the solution to K-S which 
vanishes at r = 0. 
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In a forthcoming paper we will incorporate the time derivatives of the temper- 
ature and enthalpy in the model; the front equation will be more involved and of 
higher order in time, as in [1]. Another issue we intend to address is the derivation 
of the front equation as a solvability condition in the spirit of O Hj . 

2. Some mathematical setting 

In this section we introduce some notation, the functional spaces and operators 
we will use below. We will mainly use the discrete Fourier transform with respect 
to the variable y. For this purpose, given a function / : (— £/2, £/2) C, we denote 
by f{k) its fc-th Fourier coefficient, that is, we write 

f{y)^Y.^{k)wu{y), ye {-1/2,1/2), 

k=0 

where {wk} is a complete set of (complex valued) eigenf unctions of the operator 

A - Dyy : D{A) = H^~i/2,i/2) L''{~£/2,£/2), 
with ^-periodic boundary conditions, corresponding to the non-positive eigenvalues 
47r2 47r2 IBtt^ IGtt^ SGTr^ 

Z^"' £2 ' £2 ' £2 ' ■ ■ ■ 

We shall find it convenient to label this sequence as 

- -AoW > -Xi{e) = > -A3(^) = -X4{e) >■■■ 

When there is no damage of confusion, we simply write Afc instead of Afc(^). 

When / depends also on t and/or x, by /(•, k) we denote the k-th Fourier coef- 
ficient of / with respect to y. For instance, for fixed t and x, f{t,x,k) will denote 
the fc-th Fourier coefficient of the function f{t, x, •). 

For integer or arbitrary real s, we denote by i?| the usual Sobolev spaces of 
order s consisting of ^-periodic (generalized) functions, which we will conveniently 
represent as 

{+00 -1-00 1 

w = Y^ akWk ■■ Xlal < +00 I , 
fc=0 fc=0 J 

with norm 

-t-CXD 

||«;||2 = ^A^4. 

k=0 

For fc = 0, we simply write instead of and | • I2 instead of |1 • |lo- 

We recall that for any /3 > and 7 e (0, 1) the operator (/ — /3A)'^ has H^"' as a 

domain and it is defined by its symbol ((1 -f PX^)'^) (see e.g., [TU Thm. 4.33]). 
Next, for any n = 0, 1, . . . and any /3 G [0, 1), we set 

C;+^ = {/ e C-+^{[-e/2,£/2]) : f^'H~e/2) = f^'\£/2), J < n}. 

C^'^^ is endowed with the Euclidean norm of C"+'^([— ^/2, i/2]). Finally, we denote 
by II • Hoc the sup-norm. 



8 



C.-M. BRAUNER, J. HULSHOF, L. LORENZI AND G. I. SIVASHINSKY 



(3.1) 



3. The derivation of a self-consistent equation for the front 

The aim of this section is the derivation of a self-consistent equation (in the 
Fourier variables) for the front For this purpose, we rewrite Problem (|1.13p . 
making 6 and S explicit. We get 

Vx — A(u — au) — a{ipt + {^yY){x + l)e^ — acpyyxe^, x < 0, 
- Av ^ 0, X > 0, 

uiO) = [v] - 0, 
^'(0)-u,(0) = i((^^,)^ 
. [vx] = -au.x{0). 

In what follows, we assume that {u,v,ip) is a sufficiently smooth solution to Prob- 
lem (13. ip such that the function x i-^ e~^^^u{t, x,y) is bounded in (— oo,0] and the 
function x i-^ e~^^^v{t, x,y) is bounded in M. As it has been stressed in the Intro- 
duction, we use the first equation in (13. ip and the boundary condition u{-, 0, •) = 
as a pseudo-differential Stefan condition. We solve the problem for u via discrete 
Fourier transform. This leads us to the infinitely many equations 
(3.2) 

Ux{t, x, k) - Uxx{t, X, k) + \ku{t, X, k) = [(pt{t, k) + {^PyYit, k) + \k(p{t, k)^ e^, 

for k = 0,1,2,..., where we recall that — Afe = — Afe(^) is the fc-th eigenvalue of 
the realization of the operator Dyy in L^. For notational convenience we set Vk = 
^ -I- ^Vl + 4Afc for any fc = 0, 1, . . . A straightforward computation reveals that the 
solution to (13.211 which vanishes at x 



and tends to as a; — s- — oo not slower 
than e^^/^ is given by 

u{t, X, 0) = - (^t(t, 0) + {^{t, 0)) xe\ x< 0. 

Let us now consider the problem for v, where we disregard (for the moment) the 
condition 0, •) — M2;(-, 0, •) = ^{(py)'^. Taking the Fourier transform (with respect 
to the variable y), we get the Cauchy problems 

' (t, X, 0) - [t, X, 0) = a (<^t [t, 0) + i^pyj^it, 0)) {2x + 3)e", 

Vx{t,X,0) -Vxxit,X,0) = 0, 

[«(i,-,0)] =0, 

[w,(t,-,0)] = -aM,(t,0,0) =a((^t(i,0) + (^(i,0)) , 
for fc = 0, and 

Vx{t, X, fc) - Vxx{t, X, fc) + Xkv{t, X, fc) 



a; < 0, 
X >0, 



ai^x + 2-—j (^^t{t, fc) + (^(i, fc)) + aXk 
+^ (ptit, fc) + {^{t, fc) + A,^(t, fc)) e-"'^ 



1 

A^ 



Vx{t,X, fc) - Vxx{t, X, fc) + XkV{t,X, fc) = 0, 

[v{t,-,k)]=0, 

[vxit, ■, fc)] = -aux{t,0,k) = aiy~'^ {^ipt{t,k) + {ipy)^{t,k) + Xk(p{t, fc)) 



lp{t, fc)e^ 
X <0, 
X > 0, 
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for fc > 1. 

It is easy to show that 

v{t, X, 0) = -a ((pt{t, 0) + T}Pyf{t, 0)) x{x + l)e^, a; < 0, 
v{t,x,Q) = 0, x>G. 

and 

v{t, X, k) =ci,fee"'=^ + J k) + {^{t, A;)) (x + 2)e^ + a${t, k){x + l)e^ 

+ ^ (t, k) + {^{t, k) + \k(p{t, k)) a;e'^'=^ a; < 

v(i,a;,fc) =C2,fee(^-'''=)^, a; > 0, 
where 

ci,fe =T—Tr- {'^ + ^k + ^ + — ) ^(i,A;) 

1 - 2i/fe V 1 - 2z^fc jvfc / 



a / 1 2i^fe 1 1 



1 — 2vk V -^fc -^fc Afe 1 — 2z/fc 



a /2i/fc 3 1 J^fc 



1 — 2i^k \ Afc Afe Afe 1 — 2i/fc i^k , 

Now, we are in a position to determine the equation for the front. Indeed, 
rewriting the boundary condition 

1 

2' 

in Fourier variables, and using the above results, we get to the following equations 
for the front (in the Fourier coordinates): 

^t(^,0) + i(^(t,0) =0, 

H yk\^ — ) v>[t, k) 



1 - 2i/fe \ 1 - 2i/fe i/fe J Vk 
\ l-2i^fcVAfc Afc Afcl-2i^fe I'fc/ I'fc J * ' 



Let us set Xk = \/l + 4Afc. Then, the equation for reads (in terms of Xk) as 
follows: 

^xf ^) + Xl(Xk + l) 



2X2(Xfe + 1) 
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for any fc = 0, 1, 2, . . ., or, equivalently, 



(3.3) #,(t, k) = '-0it, k) + ^ ' - {iPy)Ht, k), 

A| + aXk - a "^fe ^ ^ 



1 - - a) , , , Xi - iXl - AaXk + 4a ^ 

1 



{Xi + aXu ~ a)ipt{t, k) --(1 - Xi){Xi - a^t, k) 

+ - iXl - 4aX, + 4a) (^(t, k) 

= {-A\l + {a-l)\k)(p{t,k) 

+ \{Xl - iXl - 4aX, + 4a) (^(t, k), 
for any fc = 0, 1, . . . Hence, we can say that (p solves the equations 

(3.4) |^(^)=^(^) + ^((^,)2) 
and 

(3.5) = + m~\^{[^yf) ^(^) + ^a^yf), 
where the operators i^, 5^ and are defined through their symbols 

(3.6) hk = ^fe + aXk - a, 

(3.7) Sfe = -4A^ + (a-l)Afc, 

(3.8) /fe = i(X3-3X2-4aXfe + 4a), 
for any fc = 0, 1, . . . 

4. Stability of the front 

In this section we are interested in studying the stability and instability properties 
of the null solution to the Equations (|3.4p and p. 51) . In this respect we need to study 
the symbols appearing in p.3p . 

4.1. Study of the symbols. In this subsection, we study the main properties of 
the operators y, and whose symbols are respectively defined by 

(EH)- (EH) and by 

il-Xi)iXi~a) 

Lk ^ 



gk 



A{Xl + aXk - a) ' 

Xl ~ ZXl - 4aXfc + 4a 
4(^2 + aXk ~ a) ' 



for any fc = 0, 1, . . . Even if all these operators depend on a, we prefer not to stress 
explicitly the dependence on a to avoid cumbersome notations. 

Proposition 4.1. The following properties are met. 

(i) The operator ^ admits a realization L in which is a sectorial operator. 
Moreover, its spectrum consists of the sequence (Ik)- In particular, is a 
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simple eigenvalue of L. The spectral projection associated with this eigenvalue 
is the operator 11 defined by 



Finally, a-{L) \ {0} C (— oo,0) if and only if a < Uc {see (|1.17p ). 

(ii) The operator SS admits a bounded realization B mapping into . More- 
over, B it is invertible. 

(iii) The operator admits a bounded realization F mapping into . 

(iv) The operator admits a bounded realization G mapping into . 

(v) The realization of the operator .5^ in L'^ is the operator 

S = -ADyyyy + - l)Dyy, 

with as domain. 
Proof, (i). To begin with, we observe that 

^ Afc(4Afc + l^Q) 

aV4Afc + 1 + 4Afe + 1 - a ' 

Hence, we can split 



_ aAfcVl + 4Afc 

Ik — — Afe H — ; — Afc + ii fc, 

aVl + 4Afc + 4Afe + 1 - a 



for any k — 0, 1, . . . Note that li^k ~ §v ^fc as fc ^ +oo. Hence, from the above 
splitting of the symbol (Ik) it follows at once that the operator ^ admits a realiza- 
tion L in with domain D{L) = which can be split as L = A + Li, where Li 
is a bounded operator from into L^, and A is the realization of Dyy in with 
domain H^. Since is an intermediate space of class J1/2 between and D{A), 
[131 Prop. 2.4.1(i)] applies and shows that L is sectorial. 

Let us now compute the spectrum of the operator L. For this purpose, we observe 
that, since D{L) is compactly embedded into L^, cf{L) consists of eigenvalues only. 
Further, if A is an eigenvalue of L, then there exists a not identically vanishing 
function -0 such that = A'0. In the Fourier variables, the previous equation 
leads to the infinitely many equations 

\i}{k)~lki>{k) =0, A: = 0,1,2,... 

If A ^ Ik, then '(/'(fc) = 0. Hence, if A is not an element of the sequence [It), A is in 
the resolvent set of L. On the other hand, it is clear that the sequence {Ik) consists 
of eigenvalues of L. So (j{L) = {Ik : = 0, 1, . . .}. 

Since ^ — cx) as fc — s- +00, is an isolated point of the spectrum of L and the 
corresponding eigenspace is one-dimensional. Let us prove that H is the spectral 
projection associated with such an eigenvalue. For this purpose, we prove that 
is a simple pole of the function A 1-^ i?(A, L) and compute the residual at 0. Note 
that for any A ^ cr{L) and any ij] e it holds that 

+00 
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Hence, 

+00 ^ +00 ^ 

k=l k=l ^ 

Hence, for lAI < i min we can estimate 

^ fc=l,2,... 



iAi?(A,i)^-nv^i^<^ ^ ' ^2n...2 ^ m'' ^ ,2n...2 ^ m'' 



fc=i 



A-;*. 



^min ^min 



where ^min = niin„=i^2,... |^n| > 0. This shows that R{XI — L) has a simple pole at 
A = and its residual is the operator H, which turns out to be spectral projection 
associated with the eigenvalue 0, which is simple. For more details, we refer the 
reader to e.g., [HI Prop. A.1.2 & A.2.1]. 

To conclude the proof of point (i), we observe that ^fe < 0, for k> 1, if and only if 
1 + 4Afc — a > 0. Since (Afe) is a nondecreasing sequence, < for any A: = 1, 2, . . ., 
if and only if 4Ai + 1 — a > 0, i.e., if and only if a < ac- 

(ii), (iii) & (iv). It is enough to observe that hk ^ 4Aa.., fk ^ 2A^^^, gk ^ 



as fc — > +00 and bk for any fc = 0, 1, . . . 

(v). It is immediate and, hence, omitted. □ 

4.2. Proof of Theorem 11.11 The proof is rather classical and is based on the 
results in Propositions 14.11 Nevertheless, for the reader's convenience we go into 
details. We split the proof in two steps: in the first one we deal with Equation (|3.5p 
and in the second one we consider Equation (|3.4[) . 

Step 1. Using classical arguments based on a fixed point argument, one can show 
that for any a G K. and any T > 0, there exists ro > such that, if HvolU < ''o, the 
Cauchy problem 

^t{t,v) = {L^{t,-)){y) + {G{{vy{t,-)f)){y), t > 0, |y| < |, 

ip{0,y) ^ Lpo{y), |y| < i 

admits a unique solution ^ e Uee(o 1) ■^6{T), where 



(4.1) 



.reiT) = e C{[0,T];H^)nC\[0,T];L^) : sup £'[V']c«([.,t];H? 

L 0<£<T « 

This can be proved slightly adapting the proof of Thm. 8.1.1]. The crucial 
point is the estimate 

(4.2) s'lGii^yit, - G((^,(,s, .))')|2 < C\s'\^yi.s, ■)mc^[s,nL^)\t - s|^ 

for any < s < t < T, some positive constant Ci and any ip £ ^e{T) [9 e (0, 1)). 
To prove this estimate it suffices to observe that, by Proposition 14. iT iv) 

\G{{^y{t,-)f)-G{{Us,-)f)\l 

<C\{i^y{t, - i^Pyis, ■))% + GlDyi^Pyit, •))' - D y y {s , 

<C\l(jy{t, •) - i'vis, ■)\2\\i'v{t, •) + i^yis, ■)\\oo 

+ C\^yyit, ■)\2Hyit, •) - iJy{s, OIU 
+ CHyis, ■)\\c^\'^yyit, ') - i^yyis, Ob 
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<C\'ipy{t, •) - lPy{s, ■)\2\lljyy{t, •) + Ipyyis, ■)\2 

+ Cillbyyit, ■)\2 + C'lM-"^., ■)\2) Hyy{t, •) - ^Jyyis, Oh, 

for any < s < T, where the last side of the previous chain of inequahties follows 
from Poincare-Wirtinger inequality, and C denotes a positive constant, independent 
of s, i and -0, which may vary from line to line. Estimate (|4.2p now follows at once. 

Let us now prove properties (a) and (b). It is convenient to split the solution tp to 
Equation ([XS]) along Il{L'^) and (/-n)(i2). We get if{t, y) = p{t)wo+ip{t, y) for any 
t > {) and any y E [—£/2, £/2]. Since 11 commutes with both the time and the spatial 
derivatives, Il{Dtip) = Dtll{ip) = p' and (/ - Il){Dyip) = Dy{I - Il){ip) = DyW. 
Moreover, for any ijj G H^, G{ip) — X^iS 5fcV'(^)'^fc) so that 

nG(V) = goV^(o)--inv.. 

Hence, projecting the Cauchy problem (|4.ip along n(L^) and (/ — n)(L^), we get 
the two self-consistent equations for p and ip: 

r p'{t) = -^u{{^y{t,-)r), t>o, 

1 p(o) = n(^o). 



(4.3) 
and 



(4.4) 



Mt,y) = {Lij{t,-)){y) + {I-U){G{{i^y{t,-yf)){y), t > 0, |y| < f , 
^it,~i/2)^^it,i/2), t>0, 

4'y{t,-e/2) = ^y{t,e/2), t>o, 

Clearly, the stability of the null solution to Equation (|3.5|) depends only on the 
stability of the null solution to the equation — Lt/j + {I ~ H){G{{tpy)'^)), set in 

(/-n)(L2). 

Note that the part of the operator L in (/ — n)(L^) is still a sectorial operator, 
and its spectrum is (j{L) \ {0} = {/^ : fc = 1, 2, . . .}. In particular, all the elements 
of (7{L) \ {0} lie in (— oo, 0). Hence, the linearized stability principle applies to this 
situation. More specifically, in the case when a < ac all the eigenvalues of the 
part of L in (/ — H)(i^) are contained in the plane {A G C : ReA < 0}. Hence, 
up to replacing ro with a smaller value (if needed), for any ip g B{0,r) C i^, the 
solution ip to Problem ()4.4p exists for all the positive times. Moreover, for any 
uj > max{/fc : fc = 1, 2, . . .}, there exists a positive constant Cui such that 

\Mt, ■)\2 + ■)\\2 < C^e^^hoh, t > 0. 

As a byproduct, we can infer that the solution to Problem ()4.3p exists for all the 
positive times and 

lim p{t) = Poo ■■= Ti^o-l I n((V'y(t, ■)f)dt. 

Coming back to Problem (|4.ip . the above results show that, if a < ac, this 
problem admits a unique solution, defined for all the positive times. Moreover, 

Wt{t, Ob + \W{t, ■) -Poolloo + ||</'y(i, OIloo + Wyy{t, Ob < ^Le"* || (/Sq || 2, 

for any t > 0, any uj as above and some positive constant P^^ independent of s, (po 
and ip, i.e., the null solution to Equation (|3.5p is (orbitally) stable with asymptotic 
phase. 
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In the case when a > ac the spectrum of L|(/_n)(i2) contains (a finite number of) 
eigenvalues with positive real part. Hence, the equation -0t = L%p+{I—U){G{{^py)'^)) 
admits a backward solution, exponentially decreasing to at —oo and this implies 
that the null solution to Problem (|4.4p and, consequently, the null solution to Prob- 
lem (|4.ip are unstable. For further details, we refer the reader to e.g., [5] and [T51 
Thm. 9.1.2 & 9.1.3]. 

Step 2. We focus on the case when a < ac, the other case being simpler. Of 
course, we just need to deal with the function ijj = {I — n)(^. We assume that 
(^0 G H^- We are going to show that for any uj G (0, maxfe=i,2,... h), it holds that 

supe"'^*||(^(i, •)||4 + supe"'^*||(/?t(<, Olh < +oo. 

t>0 t>0 

For this purpose, let us consider the differentiated problem 
(4.5) 

Pt{t,y)^{Lp{t,-)){y) + {Dyy{I-n){G{{3^{p{t,-))f)){y), t > 0, \y\ < 
Py{t,~i/2)^ Py{t,e/2), t>0, 

for the unknown p = ^Jyy Here, 

(4.6) ^(C) = {I-n)(^y^ av)dv^ , C G i'. 

This problem has the same structure as Problem (|4.4p . and, by assumptions, 
Dyyipo e H^. Therefore, up to taking a smaller tq (if necessary), if ||Dj,j,(/3o||2 < ^Oi 
Problem gH) has a solution p which belongs to C^{[0,T]; L'^) D C{[Q,T];H^) for 
any T > 0. Moreover, 

supe~"'*||p(i,-)||2 < +00, 

t>0 

and p(t, •) = (/- n)p(i, •) for any t > 0. 

Let us show that ip{t,-) — ^'^{p{t,-)) for any < > 0. Clearly, the function 
^ = 3^^{p) belongs to C([0, +oo); H^) n C^HO, +oo); L^). Moreover, it belongs to 
^i/2{T) for any T > 0. Indeed, belongs to the class J1/2 between and H^. 
This means that 

w^it, •) - ^{s, ■)h < ci*(<, •) - ^{s, Oil w^it, •) - ^{s, oiii 

for any < s < t < T and some positive constant C, independent of s, t and 
From this estimate, it is clear that e C^/'^{[0,T]- Hf) C S"i/2(r) for any T > 0. 

Further, Dyy* = p and A* = ^'^{Dtp), so that DyyDt'^ = DtDyy"^ = pt. It 
turns out that 

(i) i?j,j,(A^' -£*-(/- n)G((*^)2)) = 0, {ii) *(o, •) = (/- n)(^o- 

Hence, Dt^ - L'^ -{I -\l)G{{'^ yY) = a{t)+b{t)y for some functions a, 6 : [0, +00) 
R. Since ^f*, and G((*j/)^) are continuous functions in [0, +00) x [-£o/2,4/2] 
and are periodic with respect to y, it follows that Dt^ — i^t — {I — H)G((^j^)^) 
is periodic with respect to y as well. Moreover, this latter function belongs to 
(/ — H)(i^) since vj/ does. Hence, a = b = 0, implying that 5" and V' actually 
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coincide. We have so proved that G C{[Q,+oo):H^) and Dti> G C{[0,+oo);H'j). 
Moreover, 

supe-'^*|V'tj,y(t,-)|2 +supe-'^*||V^(t,-)||H4 < +00. 

t>0 t>0 » 

To complete the proof it suffices to show that ip solves Equation (|3.4p . but this 
follows immediately observing that tp is in the domain of both the operators B 
(see Proposition 14. iT ii) ) and S, and [ipyY is in the domain of the operator F (see 
Proposition 14. 1 f iii) ) . Further, L = B^^S and G — B~^F in H^. Since ip solves the 
differential equation ipt — Lip — G{{(py)^) — 0, applying B to both the sides of the 
equation, it now follows immediately that ip solves Equation (|3.4p . 

5. Rigorous derivation of the Kuramoto-Sivashinsky equation 
In this section we are interested in proving Theorem 11.21 

5.1. Rescaling and equation for the remainder. Let (phe a solution to ()1.16|) . 

We set a = I + e and define the rcscalcd dependent and independent variables: 

The spatial period is now = £q/^, for some Iq > A-k fixed, see the Introduction. 
A straightforward computation reveals that the function satisfies the equation 

(5.1) |.^^(^)=j5-(^)+^^((7^^^)2)^ 

where 

= / - 4el?^„ + (1 + (v/^-4£A,„ - 1) , 

— —ADfif-jf-jf-j — Drjri-) 

^e = \[{.I~ 4el?^^)i - 3(/ - AeDr^r,) - 4(1 + e) ( ^1 - AeDr,^ - /) } {Dr,-f . 

Note that, if we denote by (Ajt(^)) the sequence of the eigenvalues of the second- 
order derivative with periodic boundary conditions in [— ^/2,£/2], it turns out that 
Afe(^e) = eAfc(^o) eAfc, for any fc = 0, 1, . . .. Hence, the symbols of the operators 
^e, y and are 

= Xlj^ + (1 + e)X,,k -l-e, 
Sk = -Afc(4Afc - 1), 

(5.2) Uk = liXl, - 3X1, - 4(1 + e)X,,,. + 4 + 4e), 
for any fc = 0, 1, . . ., where 

Xe,k = v/l + 4eAfc, fc = 0,l,... 
Hence, the equation for the function -ip (in Fourier coordinates) reads 

be^kMr, k) = -Afc(4Afe - 1)^{t, k) + Wv^(r, fc), 

for any fc = 0, 1, . . . Note that the leading terms (at order in e) of 6e_fe and f^^k are 
1 and —1/2, respectively. Hence, at the zero-order, we recover the K-S equation 

<I> + 4$ + <I> + "1^=0 
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As we remind it in the Introduction, this equation has been thoroughly studied by 
many authors. For our purposes, we need the foUowing classical result. For the 
reader's convenience we provide a rather detailed proof in Appendix lAl 

Theorem 5.1. Let $o G for some m > 4 and fix T > 0. Then, the Cauchy 

problem 

(5.3) 

(r, r;) = -4$^^^^(r, r/) - $^^(t, i]) - i($^(T, ^)f , r > 0, |»7| < 
D^$(t,-4/2) = i?^$(r,^o/2), T>0, fc = 0,l,2,3, 

[ $(0,77) = $0(77), I77I < f, 

admits a unique solution $ e C([0,r];i?|") such that $^ e C([0, T]; iJ"""). 
The above (heuristical) arguments suggest to split as follows: 

Tp = <^ + epe. 

To avoid cumbersome notation, we simply write p for p^, when there is no damage 
of confusion. By assumptions (see Theorem 1 1.2p . the initial condition for p is 

p(0,-)=0. 

Replacing into (|5.ip we get, after simplifying by e, 

(5.4) ^^e(p) + ^e($r) = ^(p) + .^s((*,)') + e^s((p,)') + 2^,i%p^), 

for any fc = 0, 1, . . ., where the symbols of the operators and are 

(5.5) he,k = + (1 + £)Xe,k ~2-e), 

(5.6) me,fe = ^{X^k - 3^efe - 4(1 + e)X,,k + 6 + As), 
for any fc = 0, 1, . . . 

Proposition 5.2. Fia; e e (0, 1]. Then there exists a positive constant C* such that 
the following properties are met: 

(a) for any s = 2,3,..., the operators SS^ and admit bounded realizations 
and Hf,, respectively, mapping into H^^^ . Moreover 



for any e G (0, 1] and any s as above. Finally, the operator is invertible both 
from to for any s = 2, 3, . . .; 

(b) for any s > 3, the operators ".nd admit bounded realizations Fg and Mg, 
respectively, mapping H" into H'^^^ . Moreover, 

for any e € (0, 1] and any s = 3, 4, . . . 
Proof, (a). A straightforward computation shows that 

\he,k\ = 4Afc + 4(£ + l)Afc ^ ^ + i)Afe = (6 + 2£)Afe, 
Vl + 4eAfe + 1 
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for any k = 0, 1, . . . and any e G (0, 1]. This shows that admits a bounded 
reahzation mapping into for any s > 2 and its norm can be bounded by 

a constant, independent of e e (0, 1]. 

Since b^^k — £^£,fc + 1 for any fc = 0, 1, . . ., the boundedness of the operator 
from to H^^^ follows at once. 

Showing that the operator is invertible from into H^~^ is an easy task. 
It suffices to observe that b^.k > 4:eXk + 1 for any fc = 0, 1, . . .. 

(b). Since fk = snis^k — 1/2 for any fc = 0,1,..., we can limit ourselves to 
considering the operator A simple computation shows that 

\m.M < + + 3A. + 4i±£ (v/TTiiA^ - 1 

- (l + 4eAfe)i+l 

lee^Af, 12eA2 

< ^ H ^ + 22Afe 

"(l + 4eAfe)i l + 4eAfe 

< 2^/iA| + 25Afc, 

for any fc = 0, 1, 2, . . . Hence, is well defined (and bounded) in i7| with values 
in H^~^ for any s > 3. Since its symbol can be estimated from above uniformly 
with respect to e G (0, 1], the assertion follows immediately. □ 
Since all the operators appearing in (|5.4p commute with D^, the differentiated 
problem for C := reads as follows: 

(5.7) -^^eiO + ^e(*r) = ^(C) + ^e((*'),) + e^e((C'),) + 2^e((*C),), 

where we have set ^t* = Obviously, it has a null initial condition at time r = 0. 
For simplicity, we denote by D. For an integer n > 1, is the differentiation 
operator of order n. We also set £)" = Id. 

5.2. Formal a priori estimates. For any rt = 0, 1, 2, . . . and any T > 0, we set 

^niT) = {C e Ci[0,T];H^''^^^)nC\[0,T];L') : (r e C([0, T]; i?,'^("+^))} , 

where a V 6 := max{a,6}. The main result of this subsection is contained in the 
following theorem. 

Theorem 5.3. Fix an integer n > andT > 0. Further, fixm in Theorem ] 5. 1\ large 
enough such that * e C([0, T]; i/j"+^) n C^{[0,T]; L^) and e C{[0,T]; H^+'^). 
Then, there exist ei = ei{n,T) e (0,1) and Kn = Kn{n,T) > such that, if 
C S ?Ki{Ti) is a solution on the time interval [0, Ti] of Equation (|5.7p for some 
Ti <T , then 

(5.8) sup / |i?"C(T,.)|'d77 + e sup \D-+\{T,-)\^dr^<Kr,, 

whenever < e < ei. 

To prove (|5.8p . we multiply both sides of the equation (|5.7p by (— 1)"_D^"(^ and 
integrate by parts over (— ^o/2, ^o/2). We thus get 

r— r— 

(-1)" / ]^ S,(C.(r,.))i?'"C(r,-)d^- (-1)" / I S{C{r,-))D^-ar,-)dv 
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- (-1)" J ^ (i/e(*r(r, •)) - Ms((*'),(r, •))) D'-Cir, ■)drj 
+ (-!)"£ / ^ F,((C^),(r,.))i?^"C(r,-)rfry 



(5.9) +2(-l)" / ^ F,{{-^0,{T,-))D^-C{T,-)dri, 



where S — —'^Drjnriri — D^^. 

In the foUowing lemmata we estimate all the terms appearing in the previous 
equation. We first deal with the left-hand side of (|5.9p which consists of the "benign" 
terms. 

Lemma 5.4. Fix n = 0,1, e > 0, Ti < T and ( e ^„(Ti). Then, 

(-1)" / ] i?.(C.(r,-))i^'"C(r,-)d^-(-l)" / ] 5(C(T,-))^'"C(r,-)d^ 



2dT 



«0 

d ~- 
2 dr 



(i?,I?"C(r,-),i?"C(T,-))d'7 

<0 



(5.10) + 4 r |i^"+'C(r, OPrf^ - / ' |i?"+'C(r, OP^^, 

where : — > is i/ie positive operator whose symbol is {X^,k ~ 1)- 
Proof. For any ^ g ^%(Ti), we can estimate 



/ i?e(Cr(T, •))i?'"C(r, O^^'/ = (-1)" E KkCr{T, k)XiaT, k) 

+ 00 +CXD 

^(-1)" 5^ A^C.(r, mr, k) + 4(-l)"eE AI!+iCr(T, A;)C(r, k) 
+ (-1)"(1 + e)J2 >^kiXe,k - l)Cr(r, k)l{T, k) 



fc=l 



C.(r, •)i^'"C(T, ■)dv - 4e / Cr,,(T, ■)D^^aT, ■)dr, 



+ {l + e) Cr{r, - AeD,„ - I)D^"C{r, ■)drj 

= (-!)" / (i?"C)r(T,-)I>"C(T,-)d^? + 4(-l)"e / „ D"+l(Cr(T,-))I>"+^C(T,-)d^ 



+ (-l)"(l + e) / ^^(i?ei^"C(T,-),(^"C)r(r,-))d^ 
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(5.11) 



for any r G [0, Ti]. On the other hand, a straightforward computation shows that 



(5.12) 



5(C(r,.))i?'"C(r,-)d^=-4(-l)" / , |i?"+\(r, .)N^ 



+ (-1)" / |2?"+lC(T,-)Pd^, 



for any r e [0,Ti]. Combining (|5.1ip and (|5.12p . Estimate (|5.10p follows at once. 
□ 

We now deal with the other terms in (ISj 



Lemma 5.5. Fix n = 0,1, . . ., Ti < T and assume that ^ £ C{[0,Ti];H^+'^) and 
^'t G C([0, Ti]; iJj"^^). Then, there exist a positive constant Cn, independent of 
e e (0,1] and Ti, and constants K[{n,'i') and 7^2(71, 'J) such that the following 
estimates hold 



(5.13) 



(i7e(*r(r, v)) - M,((vl/2),(r, •))) i?^"C(T, 



<K{{n,^) + \D-ar,-)\l; 



(5.14) 



FM')v)D'''ar,-)dv 



< CnS-^\D-C{r, -MD^+'ar, + C„|D"C(r, ■)\t 



C^e'\D-+\iT,-)\t + C^\D-+\iT,-)\ 



(5.15) 



F,ii^CUr,-))D''Xir,-)dv 



for any t G [0,Ti] and any C G ^^i(ri). 

Proof. Fix n = 0, 1, . . . Throughout the proof C denotes a positive constant de- 
pending on n, but being independent of t, \1/ and ^, which may vary from line to 
line. 
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Estimate (|5.13p follows immediately from Proposition 15.21 Poincare-Wirtinger 
and Cauchy-Schwarz inequalities, which allow us to estimate 



{H,{^r{r. ■)) - M,((vl/2),(r, •))) i?'"C(T, ■)d^ 



{H,{<fr{r, •)) - M,{{^\iT, •))) D^ar, ■)dr) 



<K[{n,^) + \D-ar.-)\l 

for any r e [0,ri]. 

Let us now prove Estimate (|5.14p . For this purpose, we observe that 



\fe,k\ < 2V2e^Xl + SeAfe + 2{e + 1)V^A| 



k = 0,l, 



For the convenience of the reader, we explicit the splittings we use below: 



3 2+n l+n 

1 n - 1 2+n 
- + n = \ , 

2 2 2' 



l + n = 



l+n l+n 



2 2 

n — 1 n + 1 



2 2 ' 

for n > 1. (The case n = can be handled likewise with very few slight and 
straightforward changes.) Hence, for any x G C([0, Ti]; i/f ^^") we can estimate 



(5.16) 



-^A^K,fe|x;;(r,fc)|C(r,fc)| 



A.— 



+ 00 



<2V2£i K^^\^,ir, mar, k)\+3eJ2 AI+"Ix;;(t, A:)||C(t, fc)| 



fe=0 



fe=0 



+ 00 



2Vi(l + e)5]Af"|x;;(r,fc)||C(r,fc)| 



fc=0 



3 -I- /? 



fc=0 

1 

C+OO \ 2 / + 00 



/ + 00 



/ 4-00 



\k=a 



\k=0 



/ +C)C 



2 / + CXD 



2Vi(i + e) E i^r'x;;(^' E 



\fe=0 

3 + V2 



2 / + 00 



Vfc=o / \fe=o / 

=2V2£i |i?"+\(r, •)|2|^"+'C(T, + 3e|Z?"+\(T, Obl^'^+'Cl^, Oh 
+ 2V^{l + s)\D-x{r,-)\2\D-+\{r,-)U 
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V2 



|i^"x(r,-)|2|/?"+^C(r,-)|2, 
for any r G [0,ri]. 

Now, we are in a position to prove Estimate (|5.14p . For this purpose, we observe 
that, using the Leibniz formula and the Poincare-Wirtinger inequahty, it comes: 

(5.17) |i5"+'(C(r, ■))% < C{\D''^+\{r, -MD^Cir, Oh + |^"+'C(r, 

(5.18) \D-{(:{r,-))'\2<C\D'X{r,-)\l 

for any r G [0, Ti]. Replacing Estimates (|ET7)) . in ((OB)) (with x = C^), and 

using the Cauchy-Schwarz inequality and, again, the Poincare-Wirtinger inequality, 
we get 



FAieUr,-))D''X{r,-)dv 



<Cei|i?"C(r,.)|2|i?"+'C(T, 



Ce^|i?"+^C(r,-)l2l^"+'C(r,-)|2 



+ Ce\D-Cvir, •)|2|i?"+'C(r, •)|2|i?"+'C(r, Ob + Ce|i?"+iC(r, Oil 

+ CVi(l + £)|i?"C(T, Ol'l^'^+'CCr, Ob + C|i?"C(r, Ol2|i?"+'C(r, Ob 

<Ce^D'^ar,-)\2\D"+\{r,-)\l+Ce^D^+\{T,-)\l\D"+^C{r,-)\2 
+ C£|i?"+iC(r, Ol2l^"+'C(r, Ob + Ce\D"+\{T, ■)\l\D"+^C{r, Ob 
+ CV~e{l + e)\D-aT, ■)\l\D-+'ar, Ob + C\D\{r, ■)\1\D"+'C{t, Ob 

<Cei|D«C(r, Obl^"+'C(T, OI2 + C (e'|i?"+^C(r, OI2 + e|^"+'C(r, OI2) 

+ C|i?"+2C(r, 01^ + Ce'\D^+\{T, 0|^ + Ce'\D-+\{T, 0|^ + C|i?"+2C(r, OI2 

+ C\D-C{t, 01^ + Ce|Z?"+2C(r, + C|i?"C(r, OI2 + ClD^+'Cir, OI2, 

for any r G [0, Ti] and any e G (0, 1]. Now, Estimate (|5.14p follows immediately. 

To complete the proof, let us prove Estimate (|5.15p . From (|5.16p and the esti- 
mates 

|^"(*C)(t, Ob < C|i^"C(r, Obl^"*(r, Ob, 

\D-+\^0{r, Ob < C\D-+'ar, Obl^"+'*(r, Ob, 

(which can be proved using the same argument as in the proof of (|5.17p and (|5.18p ) 
we get 



2 

<V2e^D"+\{r, ■)\^\D"+\^C){r, Ob + 3£|i?"+iC(r, Obl^"+'(*C)(r, Ob 
+ 2Vi(l + e)\D-+'ar, Obl^"(*C)(r, Ob 



+ 



^±^|n"+iC(r,Obl^"(*C)(r,Ob 



<C£i|i?"+2*(r, Obl^"+'C(r, OI2 + Ce\D-+^^iT, Obl^"+'C(r, Obl^"+'C(^, 

+ CVi|i?"*(r, Obl^"C(T, Obl^"+'C(r, Ob 

+ C|i?"*(r, Obl^"C(r, Obl^"+'C(r, Ob 
<C£i|i?"+2*(r, Obl^"+'C(r, OI2 + Ce\D-+^^iT, Obl^"+'C(r, OI2 
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+ Ce\D-+H{T,-)\2\D^+\{T,-)\l 

+ C5~^\D-^{t, ■)\1\D''C{t, + C5\D-+\{t, 

for any t e [0,Ti], any e G (0,1] and any 5 > 0. Estimate (I5.15P follows taking 
C5 = 1/4. This completes the proof. □ 
We are almost ready to write the crucial a priori estimate satisfied by C,{t, •)■ For 
this purpose, we recall that 

|z?"+V|2 < |D>|j|D"+Vl2^ ^ e h;+^. 

Applying this estimate to Z?"+^C(t, •) together with Young-inequality, yields 

(5.19) |i?"+^C(r, •)l2 < l^"C(r, 012 + Jl^"+'C(r, 01^, 

for any r g [0, Ti]. Combining Lemmata 15.41 [?751 and Estimate (|5.19p allows us to 
estimate 

\^ (l^"C(r, OI2 + ^e\D-+\{T, + (1 + e)|/R:D"C(r, 0|^) 
+ (^^ - C„e ~ eK'^{n, - C„e^|i?"C(r, OI2) |^"+'C(r, OI2 

<i^l(^, + (2 + i^^(n, VI.)) |I?«C(r, 01^ + £i^^(n, nD^^'^T: OI2 

(5.20) + C„e|i?"C(T, OI2 + C^e^\D^+^C{T, Olt 
for any r G [0, Ti]. If we set 

Ae{T) = |i?"C(r, 01^ + 4£|i?"+iC(r, 01^ + (1 + e)\^,D^X{T, r e [0, T,], 

C2-4 + 2if^(n,vp), 
C3 = 2C„, 

and assume £ small enough such that 

3 

CnS + eK'^in,^) < -, 

we can rewrite Inequality (j5.20p in the more compact form 

A',{t) + (6 - 2C£iAe(r)) |7^"+\(t, < c, + c^A^ir) + c,eiA,{T))\ 

for any r G [0,^1]. 

The following lemma allows us to estimate the function A^. 

Lemma 5.6. Let Aq, cq, c\, C2, C3, e, Tq, Ti he positive constants with Ti < Tq. 
Further, let '■ [0,Ti] — > M and A^ he positive functions of class C([0,Ti]) and 
C^{[0,Ti]), respectively, that satisfy the inequalities 

A'^{T) + {co-sA,{T))f,{T)<c,+c2A,{T)+c3e{A,{T))\ re [0,Ti], 
A,(0) = 0. 

Then, there exist ei = ei(To) € (0, 1) and a constant K = K{Tq) such that Ag{T) < 
K for any r G [0, Ti] and any e G (0, ei] . 
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Proof. The proof follows basically from [1] Lemma 3.1], which deals with the case 
when = 0. Repeating the arguments in that proof, we can easily show that 
Aeir) < 4cie==^V(3c2) for any r e [0, Ti] and any e G (0,e2(To)], where £2(ro) = 

302/(160103(6==^" - 1)). 

Let us now consider the general case when does not identically vanish in [0, Ti] . 
We fix ei(ro) < £2(7o) such that 3coC2 - 4cie=2^«eo > and e e (O,eo(r)]. Since 
^^(O) = 0, there exists a maximal interval [0, T^) where Cg — eA^ > 0. We are going 
to prove that ~Ti. For this purpose, we observe that in [OjT^) the function 
satisfies the inequality A'^ < ci + C2A^ + c^eA^. Hence, from the above result it 
follows that A^{t) < (4cie=2^°)/(3c2) for any t e [0,Te], so that cq - sA^{Te) > 0. 
This clearly implies that = Ti. □ 

We are now in position to prove Theorem 15.31 Applying Lemma 15.61 it follows 
immediately that 

sup (\D"C{T,-)\l+4e\D^^+\{T,-)\l + {l + e)\^,D^\{T,-)\l) <K,,^, 
Te[o,T] ^ ^ 

for any n 0, 1, . . ., from which ()5.10p follows at once. 

5.3. Existence and uniqueness of a solution to Equation (|5.4p vanishing at 

T = 0. In this subsection we are devoted to prove the following theorem. 

Theorem 5.7. For any T > 0, there exists eo{T) > such that that, for any 
< e < eo(T), Equation (j5.7p has a unique classical solution C. on [0,T], which 
vanishes at t = 0. 



Existence part. We prove the existence of a solution C to Equation 
vanishing at 0, by a standard Faedo-Galerkin method. Let us fix ^ e (/ — H){H^ 
and expand it into a Fourier series (see Section [2]) as follows: 



+ 00 

For TV = 1, 2, . . ., we denote by Sat = PAr((/-n)(i7|)) the projection of {I-Il){Hf) 
along the vector space spanned by the functions wi, . . . , wn . 
Let us look for a solution (n G '^n to the variational problem 

r— r— r — 

(5.21) +e r FMN)n)^dr^ + ^ / ' F,{{^CN)n)^dr^. 

for all ^ G S^r. The problem is subject to the initial condition Cn{0,-) = 0. In 
terms of Fourier series, the variational formulation (|5.9p reads as follows: 

J2 bkAN{-,mk) = skCN{-,mk) 



dr 

k=l k=l 



^ {m,,fc(vI/2),,(., A:) ~ h,^k^ri-,k)}ak) 



k=l 



(5.22) +sj2fkAc%U;knk) + 2j2fe,kmNU;k)ak). 



k=l fc=l 
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Taking ^ = wj {j — 1, . . . , N) in (|5.22p . we see that the Cn{-, fc)'s verify a system 
of N ordinary differential equations with zero initial data. Hence, there exists a 
unique solution to System ()5.22p . defined on some maximal time interval [0,Tn), 
where Tn may also depend on e. 

Next, we take ^ = Civ in (|5.2ip . The estimates of Section remain valid also 
for the function ^^r. Writing such estimates, taking as Ti any number less than T/v 
and then letting T' T/v, we thus get 

(5.23) sup I \D'XN{T,-)\^dr^ + e I \D''+\N{r,-)\''dr^<Kn, 

for any n > 1. From this estimate we infer that, whenever < e < Si{T), the 
solution of the ODE system can be extended up to T. 

We now let N +oo. For this purpose, we use Estimate (|5.23[) with n = 4. It 
leads to the following facts: 

(i) the sequence {CN)NeN is bounded in C{[Q,T]; H^), with a bound possibly 
depending on e e (0, £i(T)]; 

(ii) the sequence {{CN)T)Neti is bounded in C{[0,T]; H^) with a bound possibly 
depending on e. 

Property (i) follows immediately from (|5.23p . (Note that, if ti > 5, then the 
bound is imiform in < e < £i(r).) To prove property (ii), we observe that (|5.2ip 
may be rewritten as: 

^B,{Cn) = Pn\^S{Cn) + M,{{^\) - H,{^r) + eF.iiC^),,) + 2F,((*Ca'),) 

and we use (i) and Proposition 15.21 

Now, we can make the compactness argument work for any arbitrarily fixed 
e e (0,ei(T)]. By the Sobolev embedding theorem, the sequences (CAf)AfeN and 
((CAf)T)AfeN are bounded in C([0, T]; C^^^) and in C([0, T]; C^^^^), respectively. In 
particular, by interpolation we easily see that {Dli{(N)T-)Nefi {I = 0, ...,4) are 
bounded in C^/^([0, T] x [—£q/2, io/2]). Indeed, belongs to the class Jg/g between 
C\ and C^^'^ . Hence, we can estimate 

\\Cn{t2, •) - Cn{ti, OIIc^ < IICm(t2, •) - Civ(Tl, •)l|i|ICA^(^2, •) - Civ(Tl, oiiL^ 

" '"It 

<||^rCiv||i sup \\CN{Tr)A,,\T2-Tl\K 
re[0,T] '-'It 

for any ti,T2 £ [0, T]. Since the sequence {(n) is bounded both in C^/^([0, T]; C^) 
and in C([0, T]; Cj*/^), it is bounded in C^/^{[0,T] x [-4/2,4/2]), as weh. Arzela- 
Ascoli theorem, then implies that, up to a subsequence, C,n converges in C([0, T] ; C^) 
to a function e C([0, T]; C^^^). Similarly, D'^{{(^t)n)ni£N (l — 0,1,2) converges, 
up to a subsequence, to D'^{(t-) {I — 0,1,2). Clearly, the function ( solves the 
equation (|5.4p and vanishes at t = 0. 

Uniqueness part. Assume that Ci and (2 are two classical solutions to Equation 
(|5.7p which vanish at t = 0. Then, the function x ■= Ci ^ C2 turns out to solve the 
equation 

(5.24) W = + £^e((x(Ci + C2)),) + 2F,((vE'x),)- 
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We multiply (|5.24[) by x and integrate over [— i'o/2, ^o/2]. We get 
] B,{x{r,-))x{r,-)dv^ / ' 5(x(r, •))x(t, O^?? 

+ e Fe((x(r,-)(Ci(r,-) + C2(r,-))),)x(T,-)d'? 

' F,((*x(r,-))r,)x(r,-), 

2 

for any r £ [0,r]. All the terms but f_°l)^ F.((x(r, ^(Cil^, + C2(t, •)))r,)x(T, O^?? 
have been already estimated in Lemmata 15.41 and 15.51 Hence, we just need to 
estimate this latter term. For this purpose, we observe that (I5.16P implies that 

(5.25) 

FeiixiCi+C2)Ur,-))x{r,-)dv 

<2V2el\D\x{Ci + C2))(t, ■)\2\D\{r, Ob + 3s\D\x{Ci + C2))(r, ■)\2\Dx{t, Ob 
+ 4Vi|x(T, 0(Ci + C2)(t, Ol2|i?'x(T, OI2 

+ ^^lx(r, 0(Ci + C2)(r, 0|2|i?x(r, OI2, 
for any r G [0, T]. By the a priori estimates (j5.8|) with n = 1, wc infer that 

l^(Ci+C2)(T,0l2+e|^'(Ci+C2)(r,0l2 <2ifi, re [0,T]. 
Therefore, 

lx(r, 0(Ci + C2)(t, OI2 < Ci|x(t, Ol2|i?(Ci + C2)(r, OI2 < 2CiK\x{t, OI2, 

|i?'(x(Ci + C2))(r, OI2 < + ^2)(r, Obl^^xlr, OI2 

+ Ci|i?(Ci+C2)(r,0l2|i?MT,0l2 

< 2CiKe~'\Dx{T, OI2 + 2CiX|i?2^(r, OI2, 

for any r S [0,T] and some positive constant Ci, depending on £0 only. We can 
thus continue (|5.25p getting 

^^e((x(Ci+C2)),(r,0)x(r,0rf^ 

<CKV^\Dx{r, ■)\2\D\{t, OI2 + CKe^D\{T, + CkV^\x{t, ■)\2\Dx{r, OI2 

+ CK\Dx{r, + Ck\x{t. 0|2|^x(t, OI2 

<CKe\D\{T, + CK\Dx{r. OI2 + Ck\D\{t, 

for any r g [0,T] and some positive constant Ck, depending on K only. Hence, 
combining this estimate with (|5.10p and (|5.15p yields 

1 I / M2 , ^ \T^,J^ m2 l+£ I rS-.J- m2 , lr>2. m2 



2^lx(r, OI2 + 2£^l^x(r, OI2 + ^^^1 V^-^Cr, OI2 + M,^;^ |i?^x(r, 01 
< (i^2(0, M') + CKe + 1) |x(t, OI2 + e {Ck + K^iO, ^)) \Dx{t, 
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foranyrG [0, T], where M^^k = ^-(7x6^-/^2(0, -0)£^ • Up to replacing ei(T) with 
a smaller value £0(7"), if needed, we can assume that M^.k < for any e € (0, eQ{T)]. 
Now, Gronwall lemma applies and yields C = since C(0, •) ~ 0. 

5.4. Proof of Theorem 11.21 We now return to p and to Problem (|5.4p . This 
can be done as in the proof of Theorem 11.11 The idea is simple: we look for p 
as p(t, 77) — xi^iV) ~\~ p{'t)wq., where x h^-s zero average. More precisely, we set 
X = ^(C)i where operator is defined by (|4.6p . A simple computation shows that 

B,{x) + H,i%) - Six) - A4(($„)') - sF.iixr,)^) - 2F,{^^Xv) 

is independent of 77. Since x G (/ — n)(L^), this means that 

B,{x) + H,{%) =5(x) + (/ - n)(Ma('i>,)^)) + e{I Ii){F,{{x,f)) 

+ 2(/-n)(F,($^Xr,)). 

Let us now denote by p : [0, T] ^ M the solution to the Cauchy problem 
(5.26) 

J = ~Tl{H,{%)) + Il{M,{{%f)) + eU{F,{{xr,?)) + 2n(F,(ci>,x,)), 
p(0) = 0. 

If we now set p = p + Xi immediately see that p(0, •) — p solves equation 
(153). 

Clearly, this function is the unique solution to the Equation ()5.4p which vanishes 
at T = 0. Indeed, if pi and p2 are two of such solutions, then the functions Qi := 
DriPi and Q2 '■— £'77^2 solve Equation (j5.7p and vanish at r = 0. By the above 
results, Ci and C2 do agree. This means that (/ — 11) (pi) = (/ — 11) (^2)- But then 
also n(pi) and n(p2) agree, since, as Problem ()5.26p shows, n(pi) and n(p2) are 
uniquely determined by (/ — 11) (pi). 

To complete the proof of Theorem 1 1.21 let us check that there exists M > such 
that 

(5.27) sup \p{T,'n)\<M, 

Te[o,r] 

r,£[-fo/2,fo/2l 

uniformly in < e < eo(2^)- Applying the a priori estimates in Theorem 15.31 (here 
ri = is enough) and using (|4.6p . one can easily show that 



||(/-n)(p)||,, = ll^(C)lloo < (1 + 4)7^^. 

As far as the component of p along n(L^) is concerned (which we still denote by 
p), we observe that (see (|5.2p . ()5.5p and (|5.6p ) 

n(i/,(<i'.)) = n(M,(($,)2)) = o, 
n(F,((x,)^))--in((x,)^), 



n(F,($,x,)) = -^n($,x,), 



and we can estimate 



|n(^^a(x.(r, •))'))! <Jlx.(r)|2<i sup |C(r, •)l2 < J^o, 



|n(Fa'i>,(r,-)x,(r,-)))l <\ f_l \%{r,-)x,{r,-)\dv 
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<^l'i>.(T,-)|2|C(r,-)|2 

<i/^ sup \%{t,-)\2, 

^ Te[o,T] 
for any re [0, T]. It thus follows from that 

W)\ < e f |n(^;((x.(T, •))')) + mFe{%{T, .)x,(r, 
Jo 

^ KoT+y'KoT 

re[0,T] 

for any t G [0, T]. Estimate (I5.27|) now follows immediately. 

Finally, coming back to Problem p.l6p and setting £^ — io/y/e and = T/e'^, 
one can easily conclude that, for any e e (0,eo]5 such a problem admits a unique 
classical solution ip. Moreover, 

Mt,-) - e<^>{te^,V^ ■)\\ca-ej2A/2]) < £^M, t£ [0,T,]. 

This accomplishes the proof of Theorem 11.21 



< -KoT+^KoT sup |$^(t,-)|2, 



Appendix A. Proof of Theorem 15.11 

Showing that Problem (|5.3p admits a unique solution $ e C^([0, Tq]; L^) n 
C{[0,Tq]; H^) for some Tq > is an easy task. Indeed, the operator A : — > 
is sectorial in as it has been already remarked. By [13l Prop. 2.4.1 & 2.4.4] 
the operator S = —AA^ — A \s sectorial in with domain H"^. Classical results 
for semilinear equations associated with sectorial operators show that the Cauchy 
problem (|5.3|) admits a unique solution $ with the above regularity properties. (See 
e.g., prS, Prop. 7.1.10].) <& turns out to be a fixed point of the operator F, formally 
defined by 

(F($))(r, .) = e-^$o + r ■)fds, r > 0, 

Jo 

where {e'^'^} denotes the semigroup generated by 5*. 

Using a classical continuation argument, we can extend $ to a maximal domain 
[0,Tmax) with a function (still denoted by <&) which belongs to C^([0, Tmax); i^) n 
C([0,T„ax);i?„^). 

Let us regularize $. Suppose that $o G H^- Note that S commutes with _D^. 
Hence, 

$^(t,-) = e"^(A,$o)+ / e(^-^)^A,($^(s,-))'rfs, Te[0,T„,ax). 
Jo 

Since $ G Ci([0, T^ax); i^^) n C([0, r,„ax); and belongs to the class Jk/4 
between and 77^, we can estimate 

||$(t2, •) - <i>(Ti, OIU <l*(^2, •) - <f (ri, •)ir*ll<I'(^2, •) - <&(n, oil! 

<2||<i>r||c([0,Ti];L^)ll<f ll5([0,T,];H4)|r2 - Ti|l-*, 

for any ti,T2 G [Oi^i] and any Ti < Tmax- Therefore, by the Sobolev embedding 
theorem, we can estimate 

|D^($^(t2, •))' - D^i%iTir)f\2 <\%iT2, •)|co|*^,(r2, •) - %^{tu-)\2 

+ |*w(^2, Olool^r/ra, •) - $„(ri, Ob 
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<Cti|t2 - ri|2, 

for any ti and T2 as above. This shows that D^{^^{t2, ■))'^ belongs to C-^/^([0, Ti]; L^) 
for any Ti < T^^^. Theorem 4.3.1 of (T^ impUes that $^ G C^{[0,T^^^); L^) n 
C([0,ri„ax);-ffj')- In particular, belongs to C([0, r,„ax); -ffjj^). It follows that 
= ^riT- Iterating this argument shows that, if $0 G HI" for some m G N such 



T7] 



that m > 4, then $ e c([o, r„,ax); iJf) ^^^^ ^ C([o,r„,ax);-ffr~ )' 

The rest of the proof is devoted to show that T^ax — +00. We adapt the 
arguments in [191 Thm. 2.4]. The main step is the a priori estimate 



(A.l) 



re [0,T,„ax). 



For this purpose, we introduce the function v, defined by v{T,ri) = e~ <i>^(r, 77) 
for any (r, ?/) g [0, Tmax) x [—-^0/2,^0/2]. The smoothness of $ implies that v £ 
C"'^''*([0, r,nax) X [—^0/2,^0/2]), solves the parabolic equation 



Vt — — 3W: 



ippp-l 



2lK 



(A.2) 

and satisfies the boundary conditions D^t^\{t, — ^o/2) = dII'\{t, io/2) for any r g 
[0, T) and fc = 0, 1, 2, 3. Multiplying both the sides of (|A.2p by v{t, •), integrating on 
{—£0/2, £0/2) and observing that the integral over {—£o/2,£q/2) of (u(t, •))^w^(t, •) 
vanishes for any r g [0, Tmax), we get 

In view of the estimate 



(A.3) ^\v{t, + 3\v^r,{T, - \v,{t, + 2\v{t, = 0, 



T G [0,T;„ax)- 
r e [0,Tmax), 



Formula (|A.3[) leads us to the inequality 

^\v{r,-)\l + lHr,-)\l<0, 



re [0,T„ 



from which Estimate (|A.1[) follows at once. 

We can now complete the proof. For this purpose, let us consider the function 
v]/, defined by ^'(r, 77) = <i>(T, 77) - n(cf>(T, •)) for any r £ [0, Tmax) and any ij £ 
[—£0/2, £0/2] • Applying Poincare inequality, we get 

(A.4) |$(T,-)-n($(r,-))|2 < v/^e "#"1^^*012, re[0,T„ax). 

Let us now show that the function r 1-^ n($(T, •)) satisfies a similar estimate. 
For this purpose, we fix r S [0, Tmax) and apply the operator 11 to both the sides of 
(|5.3p . Since $ and its derivatives satisfy periodic boundary conditions. 



dr 

for any r € [0, 



n($(r,.)) = n(ci>.(r,.)) 



2£o 



Taking (jA.ip into account, we can then estimate 



Hence, 

(A.5) |n($(r))| < |n($o) 



dr 



n($(r,.)) 



T£ [0,T,„ax). 



dT< |n($o)| + 
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for any t e [0, Tmax)- Estimates (|A.4p and (|A.5|) show that <& is bounded in [0, Tmax) 
with values in . Therefore, we can apply |13i Prop. 7.2.2] with 7 = 1/2, a — 1/4, 
= H^, which implies that Tmax — +00. 
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